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We study phase shps in one-dimensional topological superconducting wires. These wires have been 
proposed as building blocks for topologically protected qubits in which the quantum information 
is distributed over the length of the device and thus is immune to local sources of decoherence. 
However, phase-slips are non-local events that can result in decoherence. Phase slips in topological 
superconductors are peculiar for the reason that they occur in multiples of An (instead of 2n in 
conventional superconductors). We re-establish this fact via a beautiful analogy to the particle 
physics concept of dynamic symmetry breaking by explicitly finding a "hidden" zero mode in the 
fermion spectrum computed in the background of a 2ti phase-slip. Armed with the understanding of 
phase-slips in topological superconductors, we propose a simple experimental setup with which the 
predictions can be tested by monitoring tunneling rate of a superconducting flux quantum through 
a topological superconducting wire. 

PACS numbers: 74.20.Mn,73.63.Nm,74.50.+r 



I. INTRODUCTION 



A quantum computer, if realized, would be able to per- 
form computational tasks with an efficiency that could 
never be reached by a classical computer. Consequently, 
great effort has been put into exploring how to realize 
such a computer. One of the main challenges in do- 
ing so lies in the high sensitivity of quantum systems 
to background noise. Storing quantum information in 
topological states of matter may provide a decoherence- 
free realization of quantum computing. In particular, as 
topological states are determined by the global proper- 
ties of the system, topological qubits are expected to be 
robust to decoherence from local perturbation^. 

We focus on a specific realization of topological mat- 
ter: topological superconducting wires. To build this 
type of wire one needs to combine the properties of three 
discrete elements: a semiconducting nanowire that pro- 
vides strong spin orbit coupling, a superconducting wire 
that provides a superconducting gap via proximity ef- 
fect, and a magnetic field that opens a Zeeman gap in 
the nano-wire spectrunP''^. Topological superconduct- 
ing wires are useful for quantum computing because a 
Majorana fermion forms at the interface between a con- 
ventional and a topological superconducting wire. By 
combining several such interfaces, it is possible to create 
a topological qubit as described in Ref. 6, Further, by 
building a network of such wires, it is possible to perform 
quantum information processing by braiding the Majo- 
rana fermions, resulting in a quantum computer with 
topologically protected quantum logic gatePfC^J. 

A possible source of decoherence in such a quantum 
computer are phase-slips in the superconducting wires. 
In a superconducting ring a phase-slip fluctuation con- 
nects states with different winding number of the super- 
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FIG. 1. Schematic of a topological superconducting qubit. 
The qubit is composed of a series of conventional supercon- 
ducting wires (labeled S) and topological superconducting 
wires (labeled T). Four Majorana fermions (labeled 71 to 74) 
are located at the interfaces. We consider phase slips at three 
types of locations (labeled Li to L3). Topological supercon- 
ductors only support 47r phase slips, which can take place at 
locations Li and L3. These An phase slips do not cause de- 
coherence of the qubit. However the central S segment (at 
location L2) can support 2ti phase slips which can cause the 
qubit to decohere. See main text for details. 



conducting phase around the ring. Phase-slip are fluc- 
tuations in which the amplitude of the order parameter 
shrinks to zero at some location along the wire, which 
results in the loss of coherence between the left and right 
sides of the wire, and the phase can slip. At the con- 
clusion of the phase slip, the order parameter amplitude 
grows, and the phase coherence is reestablished. Conse- 
quently, phase slips play an important role in determin- 
ing both the dynamics of the order parameter as well as 
in determining the quantum (and the thermodynamic) 
ground state of the wire. Phase-slips can be driven by 
either quantum or thermal fluctuation [resulting in Quan- 
tum Phase Slips (QPS) or thermally activated phase slips 
(TAPS)]. TAPS tend to dominate when the temperature 
is larger than the Josephson energy for a Josephson junc- 
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tion, T > Ej (or the corresponding energy scale for a SC 
wire). In the low temperature, T < Ej, thermal fluc- 
tuations become insufficient to overcome the barrier and 
hence QPS become the dominant process. Experimen- 
tally, both TAPS^i ^3 and QP^^SEIIl have been observed 
in thin uniform superconducting wires as well as in con- 
strictions^^ and Josephson junction^^^^^. The effect of 
phase slips on topological wires has been previously con- 
sidered in Refs. l22H25l The coherent motion of vortices in 
a 2D topological superconductor is discussed in Ref. l26l - 

In this article, we investigate the effect of quantum 
phase slips on topological superconducting wires and de- 
vices. We start by discussing the consequence of phase 
slips on a superconducting qubit shown in Fig. [l] In par- 
ticular, we note that a phase-slip of 27r, which can occurs 
in the conventional superconducting wire segment, leads 
to the decoherence of the qubit while a phase slip of 47r, 
allowed in topological wire segments, leads to no deco- 
herence. As phase slips can be an important source of 
decoherence for the topological quantum computation in 
Majorana fermion systems, it is important to study such 
processes in depth. 

We consider a simplified model where phase slips only 
occur at a weak link (Josephson junction) in a topolog- 
ical superconducting wire, and construct a semiclassical 
field theory description for phase slips at the weak link. 
Then, we recover the well known fact that although the 
fermionic spectrum is 27r periodic in the phase difference 
across the w eak link, the ground state has only 47r peri- 
odicit}'^^*^^^'^. We show this in two complimentary ap- 
proaches: (1) By integrating out the fermions, the par- 
tition function becomes explicitly 47r periodic. (2) We 
show that 27r phase slips are suppressed by making an 
analogy to the concept of symmetry breaking by a chiral 
anomaly in particle physics^^ 

Explicitly, in method (2) we view a 27r phase slip as 
an instanton event in the semiclassical description. The 
amplitude of the instanton is proportional to the determi- 
nant of the fermionic kernel evaluated along the instanton 
trajectory. Following the classic calculation of t'Hooft^^, 
we explicitly obtain the eigenvalues of the fermionic ker- 
nel. We show that the spectrum contains a "hidden" 
zero mode, that we uncover by a transformation of the 
fermionic kernel into a hermitian operator, which results 
in the suppression of 27r phase slips. Motivated by this 
result, we further discuss how the suppression of 27r phase 
slips can be observed by considering the effects of phase 
slips on topological superconductors in ring geometry 
(e.g. AC SQUIDs) as well as current biased topological 
superconducting wires. 

The manuscript is organized as follows. In Sec. |ll| we 
discuss phase slips in a qubit device composed of topo- 
logical and conventional superconducting wires. Next, we 
introduce the Kitaev model of a topological superconduc- 
tor in Sec. |III| We describe, in detail, QPS in topological 
superconducting wires and identify the hidden zero mode 
in Sec. |IV[ We discuss the detection of An phase slips 



in two types of devices made of topological supercon- 
ductors: topological superconducting rings and current 
biased wires in Sec. |V[ Finally, we make concluding re- 
marks in Sec. VI The main text is supplemented by two 
appendices, in which we derive the effective action for a 
topological superconducting wire with a weak link and 
describe the discretization of the Fermion action on the 
weak link in the presence of a phase-slip. 



II. PHASE SLIPS IN A QUBIT DEVICE 

To motivate the study of phase slips in topological su- 
perconducting wires, we review a particular implemen- 
tation of a topological qubit illustrated in Fig. [l] The 
qubit is composed of three conventional superconducting 
segments and two topological superconducting segments. 
The quantum information is stored in the four Majorana 
states labeled 71 to 74. To describe how quantum in- 
formation is stored we use the basis of complex fermions 
cl = 71+^72 and cr = 73 + ^74, for the "left" and "right" 
topological segments. We can describe the state of the 
device in terms of the occupation numbers \nL^ tir) of the 
left and right complex fermions. For states of odd parity, 
we could use |0, 1) and |1,0) to represent the two states 
of the qubit. Analogously, for states of even parity, we 
could used |0,0) and |1, 1) to represent the two states of 
the qubit. 

Consider the effect of phase slips on the qubit device 
illustrated in Fig. [l] Phase slips that can potentially 
damage the quantum information in the qubit can occur 
at three typical locations labeled Li, L2, and L3. Lo- 
cations Li and L3 lie inside topological superconducting 
segments and, as we shall show later, only support 47r 
phase slips only. On the other hand, L2 lies inside a 
conventional superconductor and thus can supports 27r 
phase slips. 

To understand how a phase slip can affect the quantum 
information stored in a qubit, we appeal to the Aharonov- 
Casher effect^^. A 27r phase slip can be thought of as tak- 
ing a vortex on a closed loop trajectory around the wire, 
with the trajectory intersecting the wire at the location 
of the phase slip. The Aharonov-Casher effect states that 
when we take a flux around a charge on a closed trajec- 
tory, the wave function builds up a phase proportional 
to the charge enclosed. In particular, when a single su- 
perconducting vortex goes once around a single electron 
charge, the sign of the wave function changes. 

Let us first consider a phase slip at location Li as de- 
picted in Fig. [1] Since Li lies inside a topological su- 
perconductor, only 47r phase slips are supported, which 
is equivalent to a vortex completely encircling the right 
segment of the topological superconducting wire twice as 
depicted in Fig. |2] The double encirclement means that 
the phase of the wave function associated with fully en- 
circled fermions is unchanged irrespective of their occu- 
pation numbers according to the Aharonov-Casher effect. 
Hence, there is no overall phase accumulation related to 
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FIG. 2. Schematic of a vortex trajectory equivalent to a An 
phase shp at location Li of Fig. [l] 



the occupation number of the cr fermion in Fig. |2j 

On the other hand, the effect of a 47r phase slip on 
the quantum state when the vortex core crosses through 
a delocalized fermion, as is the case for cl fermion in 
Fig. [2j is more delicate. To work out this scenario, we 
consider a special setting where the phase slip occurs at a 
weak link. In the limiting case of an extremely weak link 
of the topological superconducting wire, there will be a 
localized fermion = ^w,l ^ '^1w,r associated with the 
weak link. Here, ^w,l{r) ^ire Majorana fermions residing 
at the left (right) of the weak link and can be combined 
with the constituent Majoranas of the cl fermion to form 
two fermions cli and cl2 that are localized to the left and 
to the right of the weak link. By the Aharonov-Casher 
effect, the wave function of cli and cl2 fermions returns 
to its initial value following a Att phase slip. Therefore, 
the Cl fermion also returns to its initial stat^^. We shall 
give explicit arguments on how this occurs for the generic 
case in appendix [A) 

Combining the results of the previous two paragraphs, 
we conclude that a An phase-slip at Li brings the qubit 
back to its initial quantum state and does not cause de- 
coherence. In a similar manner, one can argue that a 
47r phase slip at L3 does not change the quantum state. 
The only difference is that the vortex encircles the "in- 
active" cr fermion twice for a Att phase-slip at Li, which 
accumulates a phase of 27r, while it does not encircle the 
inactive fermion cl for the phase-slip at Li, which brings 
no extra phase. 

Finally, we consider the effect of a 27r phase slip at 
position L2. Again, we let the phase to the right of 
the phase-slip core wind by 27r while the phase to the 
left remains unchanged. Here we find that states with 
Cr fermion empty remain unchanged (|1,0) ll^O), 
|0,0) |0,0)), while those with cr occupied acquire a 
minus sign (|0, 1) -|0, 1), |1, 1) -|1, 1)). Therefore, 
phase-slips at L2 decohere the qubit^^. 
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FIG. 3. Schematic of a composite structure consisting of a 
semiconducting nanowire in contact with a superconductor. 
The superconductor induces pairing in the nanowire via prox- 
imity effect. The orthogonal alignment of the spin-orbit field 
e, the magnetic field and the coordinate along the wire r 
is indicated. In implementation (a) the superconductor is a 
thin homogenous wire that is susceptible to phase slips along 
it's entire length. In implementation (b) the superconduc- 
tor is rigid everywhere except a weak link, a point at which 
phase-slips can occur. 



III. SETTING: TOPOLOGICAL 
SUPERCONDUCTING WIRES 



To make concrete arguments about phase slips in topo- 
logical superconducting wires, and devices containing 
topological topological wires, we focus on the implemen- 
tation of topological superconducting wires described in 
Ref. 3, In this implementation, topological supercon- 
ductivity is not obtained as an intrinsic property of a 
material, but rather by combining various materials to 
engineer the desired properties. The main part of the 
proposed composite is a single channel semiconducting 
nanowire with strong spin orbit coupling. By applying 
a strong magnetic field, the electrons in the nanowire 
form two, well separated, spin polarized bands. Due to 
the presence of both a magnetic field and the spin orbit 
scattering, the spin polarization in the two bands is mo- 
mentum dependent. Finally, by proximity coupling the 
semiconducting nanowire to a conventional s-wave super- 
conductor, we induce p-wave pairing in the bottom band 
of the nanowire. Thus the semiconducting nanowire is 
predicted to exhibit topological superconductivity. 

In the composite implementation of topological super- 
conductivity, phase slips in the topological superconduc- 
tor are associated with phase slips in the proximity giv- 
ing superconductor. We therefore assume that the su- 
perconductor is sufficiently weak so that it can support 
phase-slip fluctuations. This can occur if the supercon- 
ductor is a sufficiently narrow wii^^^, or if there is a 
weak spot or break in the superconductor which results 
in the formation of a Josephson junction. These possi- 
bilities are schematically illustrated in Fig. [3] To model 
the composite structure, we use the Kitaev model^^J to 
describe the electrons in the semiconducting nanowire, 
and supplement it with a phenomenological model that 
describes the order parameter in the proximity giving su- 
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perconductor. 

The Kitaev model is specified by the Hamiltonian 

N N-1 



i=i 1=1 

(1) 

where N is the number of lattice sites, c] (q) is the elec- 
tron creation (annihilation) operator at site i, /i^ is the 
chemical potential at site z, t > is the hopping matrix 
element, and A^^^+i is the complex order parameter, de- 
fined on the link between sites This model can be 
thought of as the large magnetic field regime of the model 
described in Ref. HI The model supports both topological 
and conventional phases by tuning of the chemical po- 
tential, with the phase transitions occurring at = 2t. 
Thus, we can model both topological and conventional 
segments by varying fii as a function of position along 
the wire. 

To describe the dynamics of the order parameter in 
the superconductor, we need to choose whether we are 
describing a Josephson junction or a continuous thin su- 
perconducting wire. As we are interested in the effect 
of the electron degrees of freedom on phase slips, these 
details will not be especially important. In the next sec- 
tion, we shall focus on the technically simpler problem of 
phase slips at a Josephson junction (weak link). 



IV. PHASE SLIPS AT A WEAK LINK: 
HIDDEN ZERO MODE 



THE 



In this section, we construct a theory of phase slips in 
the weak link geometry illustrated in Fig. [sj^b): a semi- 
conducting wire on top of a superconducting wire with 
a single weak link. We start with this geometry as it in- 
volves fewer degrees of freedom than the continuous wire 
geometry illustrated in Fig. [sja). 

We explicitly construct an effective, low energy, model 
of the weak link geometry starting from the Kitaev 
model ([T]) in appendix [A[ From the point of view of su- 
perconductivity, the weak link geometry is a Josephson 
junction, that can be characterized by the phase differ- 
ence (/) across the weak link. From the point of view of the 
electrons in the semiconducting nano-wire, the weak link 
is a topological-conventional-topological junction. Asso- 
ciated with each topological-conventional interface, there 
is a Majorana fermion. By assumption, the weak link is 
short compared to the Fermi- wavelength in the nanowire, 
and therefore the two Majorana fermions interact to form 
a single complex fermion c^^; that is localized on the weak 
link. The low frequency effective action involves (j) and 
Cyj degrees of freedom associated with the weak link and 
is given by 



5j 



dt 



l^(dt4>f-Ej{l-cos{4>)) (2) 
+ cl {idt - Em cos((/)/2))c„ 
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FIG. 4. (a) Instanton trajectory in the sine-Gordon model, 
Eq. Q. (b) Schematic representation of the dilute instanton 
gas composed of 27r phase slips and —2n anti-phase-slips. 



In this model, the first term is phenomenological in origin 
and describes the charging energy Ec = e^/2C due to 
the capacitance C associated with the weak link. The 
Ej term describes the 27r periodic part of the potential 
energy and is primarily related to the electronic states 
of the semiconducting nanowire outside the gap. There 
can be a secondary contribution to the Ej term from 
the Josephson energy associated with the weak link in 
the underlying superconductor. The final term describes 
the sub-gap fermion c^^, localized at the weak link. The 
energy scale Em and Ej can be obtained from the Kitaev 
model, see appendix [Aj 

We begin by sketching the semi-classical dynamics of 
the phase only (i.e. sine-Gordon) model, without the 
fermionic term, as described by the real time action 



/ 



dt 



1 1 



2 8Ec 



{dt^f -Ej{l- cos(0)) 



(3) 



The potential energy associated with the second term 
of this action is 27r periodic, thus we would expect that 
the phase would be localized near 0, ±27r, ±47r .... How- 
ever, quantum fluctuations driven by the first term can 
connect these minima via phase slips. Following the in- 
stanton prescription, we can obtain a semiclassical ap- 
proximation for the tunneling matrix elemeni(^^^. The 
prescription states that we must first go to the imaginary 
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time (Euclidean) description via t ^ ir 



dr 



1 1 



2 8Ec 



(4) 



Going to the Euclidean description results in the change 
of the sign of the potential energy term. Thus, the 
minima at and 27r in the real time description, be- 
come maxima in the Euclidean description. Moreover, 
in the Euclidean description there is a classical trajec- 
tory ^ci('?") that connects these maxima: (j)ci(—oo) = 
and (j)c\{oo) = 27r, which is illustrated in Fig. Hi a). The 
instanton trajectory leads to the value of the tunneling 
matrix element, which at lowest order is 



AHt 



|27r) 



^ e 



-S(j)[(f)cl] 



(5) 



where ^'^[^ci], is the value of the action associated with 
the classical trajectory (j)c\{r). 

A complimentary approach to studying dynamics is 
to study the thermodynamical ground state. Instanton 
trajectories extremize the action and are therefore im- 
portant in the description of the thermodynamic ground 
state. Indeed, we can think of the low temperature 
ground state, associated with S^^, as a dilute gas of phase 
slips and anti-phase-slip j^^ * ^^ , which is schematically il- 
lustrated in Fig. Qb). 

At this point we are ready to ask the question of what 
is the effect of Fermions, i.e. the third term in Eq. |2j 
on the phase slips and therefore on the ground state. To 
answer this question, we consider the partition function 
corresponding to the thermodynamic ground state 



(6) 



where Sj is the Euclidean action associated with S'j . We 
are particularly interested in the low temperature regime 
T ^ 0, in which the integral in runs over a long stretch 
of imaginary time from r = 0tor = /3 = l/T. We will 
answer the question about the role of the fermions in two 
ways. First, we will integrate out the fermions and obtain 
an effective phase-only partition function that takes into 
account the contribution of the fermions. Second, we will 
appeal to a beautiful analogy to a problem in particle 
physics to show how the fermionic term breaks 27r phase 
rotation symmetry in the ground state. 



A. Method 1: Integrating out fermion 

In this subsection, our goal is to integrate over the 
fermionic degrees of freedom in the partition function 
and convert the action Eq. Q to an effective action de- 
pending only on the phase, (j). Since the fermionic part 
of the Lagrangian is quadratic, we can integrate over the 
fermionic degrees of freedom in Eq. (|6| for an arbitrary 
trajectory (j){r) and obtain the expression 



Z (X y P(^det[/C/((^)]e" 



(7) 



Here, we use the proportionality sign to accommodate the 
normalization of the fermion path integral, and /C/(^) is 
the Lagrangian density of the fermionic part of the action 



I 



drcl [dr + Em cos(^/2)] 



(8) 



subject to an anti-periodic boundary condition c^(/3) = 
-c^(0). 

To compute the fermionic determinant, we make use 
of the fact that det[/C/ ((/))] = fln'^^' where A^'s come 
from the eigenvalue problem 



JCf{(j))Un{r) = XnUn{r). 



(9) 



Solving the eigenvalue problem, for arbitrary (/)(r), we 
find the implicit expression for the eigenfunctions Un 

Un{r) = e/o"(^--^^^°«[^(^')/2])^^'. (10) 

With the anti-periodic boundary conditions, we obtain 



i7r{2n + 1) h 



(11) 



where Ii = Jq dr Em cos((/)/2) and n is an integer. Using 
a few well known identities as in Ref. 37} we now find 



det[/C/(0)]=n 

n 

= n 



i7r{2n + 1) h 
Z7r(2n + 1) 



(12) 



cosh(/i/2). 



Thus we find that the partition function becomes 
Z cx Zeff = J cosh(/i/2)e"'^^ 



(13) 



We interpret this partition function as follows. The 
fermion can be in one of two states (either even or odd 
parity), since there are no terms in the Hamiltonian 
that connect these states, the partition function splits 
into two parts: one part for even parity and the other 
part for odd parity, manifested in cosh(/i/2)e~*^'^ = 



-h/2 , 



The even and odd parity 

states are separated by the energy Em cos(0/2), and the 
effective action becomes 



-eff 



-f 

Jo 



dr 



1 1 



2 8Ec 



(dt^Y - Ej{l - cos((/.)) 
±^cos(0/2) 



(14) 



where the sign of the last term is determined by the parity 
of the fermionic state. We note that this action is called 
the double sine-Gordon model. 
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B. Method 2: the hidden zero mode 

As we have argued, the quantum (as well as the low 
temperature thermodynamic) ground state of is com- 
posed of a superposition of quantum states where (j) lo- 
calized at multiples of 27r. Due to quantum fluctuations, 
states with different 0's are connected by instantons. In 
this subsection, we explicitly show that this picture is 
significantly modified in the presence of the fermion de- 
gree of freedom, by considering the fermionic path in- 
tegral in the background of a phase slip. Indeed, what 
we find is that 27r phase slips are strongly suppressed by 
the appearance of a "hidden" zero mode in the fermionic 
determinant. As a result, the 2tt periodic symmetry of 
the spectrum is broken down to 47r periodic symmetry 
in accord to the effective action that was obtained in the 
previous section. This mechanism of symmetry breaking 
was first studied in the context of high energy physics, 
specifically it was used by t'Hooft to explain the "missing 
meson" problem of quantum chromodynamics in Ref. [32l 
see also Refs. 'SS' and 'ST, 

Consider a bounce (phase-slip followed by an anti- 
phase-slip), such that 0(0) = 27r and (j){j3) = 0. To be 
concrete, we will focus on phase slips with the functional 

form cos(0(r)/2) = tanh . In describing the rare 

instanton gas, the instantons must be separated by long 
stretches of imaginary time. Therefore, to understand a 
single instanton, we must look towards the limit /3 ^ oo. 
What do we expect in this regime? Following the above 
discussion of the partition function, we expect that the 
matrix element must be (see Ref. [33] and [34l for details) 



iHt 



|27r) 



det[/C/,27r] 
det [ICffi] 



(15) 



where /C/,27r (^/,o) is the Lagrangian density operator 
in the presence (absence) of a bounce. /C/,o is necessary 
for normalization. In what follows, we will use the simi- 
lar subscripts 27r(o) to indicate operators in the presence 
(absence) of a bounce. Specifically using Eq. 12 , the ratio 
of fermion determinants is 



det[/C/,27r] 
det [/C/,o] 



cosh 



{Uo^ cosm/2) dr) 
cosh(/3/2) 



(16) 



which becomes - e-^/^ in the limit /3 ^ oo, since with 
the bounce, the integrand in the numerator will be neg- 
ative for a large part of the interval [0,/3], and thus 

J^^cos{cP[T]/2)dT^l3. 

To uncover the "hidden" zero mode in the fermionic de- 
terminant, we first rewrite the fermionic action, Eq. ([8|, 
in a doubled form 



Sf{(l)) = j dT^^Cf{(^)^lj 





dr - Em cos(0/2) 



(17) 



V', 



where i/j^ = (c|^,c^), subjected to anti-periodic bound- 
ary conditions '0(/3) = — '^(O). Evidently, we have 
det [/C/ ((/))] = ydetjZTMI. which can be shown explic- 
itly by using the fact det[i2j] = fli where are eigen- 
values of the differential equations 



Ui{r) 



A, 



Ui{r) 



(18) 



The eigenvalues A^ can be obtained in the similar way 
as Eqs. ([9|, ([lO| and (11), and take the form A^ 



± -k for all integer n. Here, A+ correspond to 



i7r(2n+l) 

/3 _ 

Ui sector while A~ correspond to vi sector. As expected, 
the product of all A^ gives det[/C/(0)]^. 

To facilitate the analysis, we transform the differential 
operator in Eq. (18) into a difference operator Lj. By 
discretizing the interval r G [0, /3] with N lattice points, 
we first arrange the amplitudes of the wave function at 
each lattice site, Un and Vn with n G 1, . . . , A^, in a vector 
form 



(19) 



S = (ux,U2, . . . ,iXAr, Vl, ^2, . . . ,^iv)^, 



Then, the difference equation corresponding to Eq. ( |T8| ) 
becomes LjS = AS, where the difference operator takes 
the form L/ = L^^ © L^. We then have 



TU 



TV 



1 

25 
1 

26 



(20) 



where ij G 1, . . . , A/", = cos{^{n5)/2) and 5 = P/N 
is the step in imaginary time. Now, the determinant of 
the difference operator det[L/] is simply the product of 
all eigenvalues of A. 



However, discretization scheme in Eq. (20) suffers from 



the notorious fermion doubling problem and effectively 
doubles the number of fermions both for u{r) and v{r) 
sectors^^. Hence, the continuum limit of the determinant 
det[Lj] Itv^oo is not associated with det[>Cj] directly. In- 
stead, one expects the relation det[I/j] Iat^oo ^ det[i2f]^. 
By introducing the proper normalization as in Eq. mL 
we find 



det[L/,27r] 



det[L/,o] 



det[/:/,27r]^ _ det[/C/,27r] 



det[/:/,o]2 det[/C/,o] 



(21) 



We compute the spectrum of the difference operator L f 
using anti-periodic boundary conditions with constant 
(j){r) and with a 27r phase slip followed by a 27r anti- 
phase-slip, see Fig. [sja). We have to use a phase-slip fol- 
lowed by an anti-phase-slip in order to make the bound- 
ary conditions on the fermions make sense. Without 
phase-slips, the eigenspectrum of Lf^o contains two lines 
of eigenvalues in the complex plane with ReA^ = ±Em, 
see Fig. jsjb). In the presence of the phase slips, the 
eigenspectrum deforms as plotted in Fig.jsjb). However, 
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FIG. 5. (a) cos[(/)(r)/2] as a function of r for the no phase- 
slip case (blue), and a phase-slip followed by an anti-phase-slip 
trajectory (red), using /3 = 16. (b) Eigenspectrum of iv/,27r 
with antiperiodic boundary conditions, = 16, rir = 128. 
Blue dots represent the spectrum with no phase slips and red 
dots represent the spectrum with a phase-slip followed by an 
anti-phase slip, (c) Eigenspectrum of T • L/,27r, no phase slip 
on the left and phase-slip followed by an anti-phase slip on 
the right. The fermionic spectrum on the right contains four 
zero modes. 
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FIG. 6. Ratio of determinants for the phase profile pic- 
tured in Fig. [5]^ a) computed using different methods. (1) 
following the prescription of Method 1 we integrate out the 
fermionic degrees of freedom without discretization, and raise 
the final answer to the fourth power to compensate for the 
two fermion doublings in the discretized methods (labeled: 
cosh(/3/2)~^). (2) following prescription of Method 2, we 
compute the fermion determinants on a discrete lattice (la- 
beled: |L/,27r|/|iv/,o|). (3) following Method 2, by con- 
structing the ratio of the four smallest fermion eigenvalues 
(Ai,27r/Ai,o) X ••• X (A4,27r/A4,o) of i^/,27r and respec- 
tively (labeled: zero modes). Comparison of the three curves 
indicates that the suppression of tunneling is indeed controlled 
by the zero modes, with the small offset being a non-universal 
feature associated with the duration of the phase slip. 



in the presence of phase-slips, the spectrum contains no 
obvious zero modes. 

The final step needed to uncover the zero mode is to 
consider the operator Hf = T ■ Lf, where T = ia^ ^ In 
and Iat is a X identity matrix. We note that 
this transformation does not change the determinant, 
det[i^/] = det[Lj] (up to a sign, which gets cancelled in 
the normalization). While the operator Lf is not hermi- 
tian, the transformed operator Hf is hermitian. Indeed, 
the eigenspectrum of the Hf operator without phase slips 
looks like a gapped spectrum, with the gap set by Em^ 
see Fig. [sjc). On the other hand, for the phase-slip fol- 
lowed by an anti-phase-slip <p{r) trajectory depicted in 
Fig. [sja), we find that the gap is occupied by four modes 
with near zero eigenvalues. As the splitting of these 
modes from zero depends exponentially on the separa- 
tion of the two phase-slips, we shall refer to these modes 
as the zero modes. 

On closer inspection, the Hf Hamiltonian looks like the 
Hamiltonian of polyacetylene. In continuum notation, 
the operator Hf is 

rr _ f -S,-EmCOs((/)/2)\ 

[dr- Em cos{^/2) ; ^^^^ 

where r represents the position along the polyacetylene 
chain. Now, we can leverage the well known properties of 
the polyacetylene Hamiltonian to understand our Joseph- 
son junctions action: each time the mass changes sign 
(i.e. <j) phase slips by 27r) there appears an extra zero 



mode that is localized on the kink (phase-slip). Because 
of Fermion doubling, in the discrete version we actually 
find two zero-modes associated with each kink. In case 
there is more than one kink, the zero modes will be split, 
with the splitting being exponentially suppressed in the 
separation of the kinks. Indeed, in Fig. |5|c) we see a 
signature of this effect, with four zero modes appearing 
in the gap, once we introduce two kinks (a phase-slip fol- 
lowed by an anti-phase-slip). In summary, going back to 
the original undoubled model Eq. ([2|, each phase-slip is 
associated with 1/2 zero mode. 

We pause to remark on the relation between the bound- 
ary conditions and the zero modes. In principle, we can 
choose open, periodic, anti-periodic or some other form 
of boundary conditions. Despite the choice of bound- 
ary conditions, each 27r phase slip will result in the ap- 
pearance of two additional zero modes in the discretized 
model. We note that for the case of anti-periodic (or 
periodic) boundary conditions, in order for the sign of 
cos(^(r)/2) to match across the boundary, phase 
slips must be added in multiples of 47r. Finally, we add 
that in order to obtain the correct value of the partition 
function, we must indeed use anti-periodic boundary con- 
ditions, see appendix of Ref. 37, 

Having found that phase-slips in the order parameter 
are associated with zero-modes in the fermion determi- 
nant, we now demonstrate that these zero modes indeed 
control the value of the fermion determinant. To test 
this, we consider the two trajectory depicted in Fig.jsja). 



First, as a consistency check, we compute the ratio of de- 
terminants for this pair of trajectories using both the 
continuum method described in the previous subsection 
and the discrete method described in this subsection. To 
make a direct comparison, we square the continuum re- 
sult in order to match the effects of fermion doubhng. 
We plot the comparison, as a function of the number of 
discretization steps in Fig. [6j The figure demonstrates 
that the two ways to compute the ratio of the fermion 
determinants converges as the number of discretization 
steps increases. Next, we compare the ratio of the de- 
terminants to the ratio of the four smallest eigenvalues, 
i.e. the product of four eigenvalues of near zero modes 
divided by the quartic of the gap, (£^m/2)^. We see that 
the ratio of the eigenvalues follows closely the ratio of 
the determinants computed using the discrete method, 
except for a small offset of order unity, see Fig. [6) The 
offset is associated with the imaginary time size of the 
phase slip. Thus the ratio of determinants is indeed con- 
trolled by the zero modes. 

In appendix |Bj we shall describe an alternative dis- 
cretization scheme for avoiding the fermion doubling 
with the cost that the spectrum of the difference op- 
erator under such scheme would not match Eq. (11). 
However, with such discretization scheme, the con- 
tinuum limit of determinant det[Lj] Iat^oo corresponds 
to det[>C/] directly. Hence, the ratio of determinant 

det[L/,27r]/det[L/,o]|iv^oo = det[/:/,27r]/det[/:/,o]. More- 
over, when diagonalizing the transformed operator Hf = 
T ■ Lf^ only two near zero modes appear around the in- 
stanton and anti-instanton in the phase field trajectory 
shown in Fig. [sja), a clear signature of the absence of 
fermion doubling. Hence, the ratio of determinants fol- 
lowing two different trajectories in Fig. (sja) is predom- 
inated by the ratio of these two smallest eigenvalues to 
the square of the gap, (£^m/2)^. 

In summary, we find that associated with a 27r phase- 
slip, there is a hidden fermionic zero mode. We can reveal 
this zero mode by transforming the Lagrangian density 
operator Lf with icr^ to find a hermitian operator Hf. 
The appearance of the zero mode suppresses 27r phase- 
slips. 



V. TOPOLOGICAL SUPERCONDUCTING 
DEVICES 

In this section we consider two different setups, which 
could be built to detect the suppression of 27r phase slips 
experimentally. One setup, shown in Fig.[7|^a), consists a 
superconductor ring interrupted by the Josephson junc- 
tion while the second setup, shown in Fig. [7|^b), is a 
normal Josephson junction with a constant supercurrent 
passed through it. While the first setup is conceptually 
cleaner as the tunneling of a flux quantum out of the 
loop is measured, the second has the threefold advantage 
that it does not involve building a loop, that it does not 
involve an inductance of a magnitude which is challeng- 



(a) 



(b) 




FIG. 7. A topological superconducting wire is place across a 
Josephson junction. In panel (a), the junction is connect to a 
superconducting ring and a magnetic flux $ can be threaded 
through the ring to bias the conductance energy. In panel (b) , 
the junction is currently biased to form a washboard potential 
that drives phase slips. 



ing to realize, and that it does not involve changing the 
inductance El but rather the bias current Is when deter- 
mining the power-law suppression of the phase-slip rate 
due to the zero-mode, see below. 



A. Ring geometry 

In the absence of the topological superconductor wire, 
the Euclidean action of the Josephson junction read^^ 



f 

Jo 



dT 



+El { m - 27r|- 



(23) 



where (j) is the gauge invariant phase difference across 
the Josephson junction and ^/^o is the ratio of the ex- 
ternal magnetic flux threaded through the ring and the 
superconducting flux quantum = h/2e. As a super- 
conductor ring interrupted by a Josephson junction is 
characterized by its critical current 1^^ its capacitance C 
and the self-inductance L of the ring, we have the fol- 
lowing energy scales: the charging energy, Ec = e^/2C, 
the Josephson energy Ej = <I>o/c/27r and the inductive 
energy El = ^q/Stt'^L. 

The potential energy is given by the last two terms of 
the action Eq. ([23|. In the absence of the inductance 
energy as in Eq. Q, the cosine potential favors states 
with (j) = 27rZ. The inductance energy breaks such de- 
generacy by favoring states with cj) ^ 27r<I>/<I>o. To still 
have well defined potential minima at ^ ~ 27rZ, we will 
assume that Ej is the largest energy scale of the sys- 
tem and hence Ej ^ El- When <I> = 0, there are a 
global minimum at (/> = and well defined local minima 
at ^ ~ ±27r. As we are interested in the occurrence of 
phase slips of 27r, i.e., tunneling or relaxation of the phase 
from one minimum to another, we can first prepare the 
system with $ = <I>o at t < such that a flux quantum 



9 



V/Ej 

1 It 


El/Ej = 0.002,- 
y 

1 \ 3 

1 . If'/'. 




S 2k 



\(b) "/"^ 


= 0.02 


\ ^ 


/ 














-2jt ( 


3 2jt 



FIG. 8. The potential profiles of V^{(l)) in the action ([25) 
are plotted in solid (red) and dotted (blue) lines for for V± 
respectively, with Em/Ej = 0.25. The dashed (green) line is 
the potential without the Majorana fermions, i.e., Em — 0. 
The panel (a) shows the typical situation for El < Em/^tv'^ 
with El/Ej = 0.002, where two degenerate minima sit at 
(/) ~ di27r. The panel (b) shows the typical situation for El > 
Em/^tt"^ with El/Ej — 0.02, where the potential minimum 
is at = and two local minima are around 6 ~ it27r. 



is trapped inside the ring and (j) = 27r. Then, we turn off 
the external flux at t = and observe the relaxation of 
phase from = 27r to which manisfests itself as voltage 
spike across the Josephson junction. 

As shown in Sec. |IV[ the low energy fermionic degrees 
of freedom of the topological superconducting wire couple 
to the gauge invariant phase difference. The effective 
action is given by 

S^= r dr^{ryl[tdr^ Em cos <l,{r)a']iP{r). (24) 
Jo ^ 

The presence of fermions influences the tunneling rate be- 
tween different phase minima. As we showed in Sec. |IV[ 
the effect of the low energy fermion can be investigated 
by two routes as detailed below. 

1. Integrating out fermions 



Following procedures in Sec. |IVA[ w e can first inte- 
grate out the fermionic action Eq. (24) and obtain the 
effective actions for <l> 



^eff 



./3 

/ dr 

Jo 



1 1 



2 8Ec 







{drCl){T)y ^Ej{l- cos ^{t)) 



+^^(/>2(r)±^cos((/>(r)/2) 



(25) 



We observe that integrating out of fermionic degrees of 
freedom simply adds the term ±Em cos(^/2)/2 into the 
original bosonic action with the choice of ± sign depend- 
ing on the fermion parity of the system. 

To understand the effective actions, we first plot the 
profiles of the potential term 



Ej{l- cos ^{r))^EL^\r)± 



Em 



cos((/)(r)/2), 
(26) 



in Fig. [s] with a shift to make all V^{0) = 0. The ini- 
tial condition is prepared such that the superconducting 
wire is at its ground state for cj) = 27r. Therefore, with 
Em > 0, the effective action should take the sector S*^, 
which will be assumed throughout the following discus- 
sions. We note that the effective potential V~^{(j)) be- 
haves qualitatively different depending on E^ is greater 
or smaller than Em/^tt'^- When the inductance energy 
is dominates. El > Em/^tt'^^ the potential has a global 
minimum at ^ = and two local minima at = ±27r. 
In contrast, when the Majorana fermion energy becomes 
substantial. El < Em/^tt'^, there are two degenerate 
minima at ±27r and a local minimum at ^ = 0. 

From the potential profiles in the El < Em/^tt'^ 
regime, we find that a phase slip from (/) = 27r to = is 
energetically unfavorable as V^^(O) > V~^{27t). Instead, 
a phase slip of 47r, tunneling between (j) = ±27r, would 
lead to a stable state. As discussed earlier, such a phase 
slip would not change the states of a qubit based on this 
system. 

For the regime where El > Em/^tt'^, an initial state 
at (/) = 27r can relax to (j) = state since now V~^{0) < 
V+(27r). The relaxation rate is given by r27r^o = Ke~^^ 
where K corresponds to the attempt rate for the tun- 
neling and Sq is the adjusted action evaluated along the 
bouncing trajectory that starts from the initial energy 
minimum ^ 27r to the bouncing point and then 
back to (pi.^^ Here, the adjusted action is defined by 
S' = — / drV~^{(j)i) such that the corresponding po- 
tential V'{(l)) = V+(^) — y+(0i) vanishes at the potential 
minimum As a rough first approximation, we can as- 
sume that K is not affected by the presence of Majorana 
fermions and plays no role for our discussion. 

To compare the relaxation rates, r2^^o (Majorana 
fermions present) and T2j^o (Majorana fermions ab- 
sent), we shall now compute the Sq for both cases. As the 
bouncing trajectory is a stationary path of the equation 
of motion, one can show that 




d4>y^V'{4>). 



(27) 



In the case of El = Em = 0, we hav e — 27r and = 0, 
and the action is S'^ = A^2Ej/Ec- When El/Ej <C 1, 
we still h ave (pj ^ 27r and ^5 ^ 0, and we can approximate 
Sq ~ A^2Ej/Ec' Qualitatively, the presence of a small 
inductance energy El/Ej <C 1 increases the relaxation 
rate only slightly, i.e., decreasing the action such that 

Sq ^ Sq\el=0' 

We observe that the suppression of tunneling rate due 
to the Majorana fermions is given by e~^^^^ where 



5Sn — So So 



0\E^ 



(28) 



is the difference between the actions. From Eq. ( |27| ), one 
can see that SS'q is of the form SS'q = ^/|^/ fy) • 
In the limit Ej ^ El ^ £^M/(47r^), one can approxi- 
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FIG. 9. The ^5*0 in Eq. (28) is evaluated numerically and 
shown in blue curve as a function of El/Ej with Ec/Ej = 1 
and Em/Ej — 0.05. The red curve is the approximate result 
shown in Eq. ( 30 ) . 



FIG. 10. The time interval Th as a function of El/ Ej is 
evaluated numerically with Ec/Ej — 1, c.f. Eq. (33), and 



shown in the blu e curve. The red curve is the approximated 
result in Eq. (35). 
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M 



.\n{Ej/EL). 



(29) 



(30) 



It is however straightforward to evaluate 5Sq numeri- 
cally, which is shown in Fig. |9] as a function of El/Ej 
with the parameter Em/Ej = 0.05 and Eq/Ej = 1. 
The red line shows the approximation result in Eq. (|30|. 
Here, the positive sign of 5Sq indicates the suppression of 
relaxation rate. In general, a smaller El/Ej and larger 
Em/ Ej leads to a stronger suppresion. We also note that 
the approximated form of / only provides a qualitative 
trend of /(§y, ^j)- However, in the following subsec- 



tion we will show that the approximate form Eq. (29) is 
indeed the fingerprint of the zero mode physics. 

2. Relation to zero modes 

In the limit that Ej ^ El ^ £^M/(47r^), we can first 
neglect the presence of the Major ana fermion and fol- 
low the bouncing trajectory of action Eq. (23). Then, 



the Majorana fermion can be integrated out with the as- 
sumption that (/)(r) follows the bouncing trajectory. Such 
a trajectory can be evaluated by realizing that 



1 



i2 



IQEc 



Ej{l 



Er 



E^ 



(31) 



is conserved along the classical trajectory. From the ini- 
tial condition, (/) = 27r and = 0, we have E = —Att^El 
and hence the classical trajectory satisfies 



d(j) 
'dt 



4y^E,{Ej{l - cos(/)) + El^^ - ^tt^El) 



(32) 



As discussed in Sec. |IVBl zero modes appear when the 
superconductor phase difference (j){r) passes through tt. 



i.e., from ^ > tt to ^ < tt or vise versa. For a bounc- 
ing event, similar to the phase trajectory depicted in 
Fig.jsja), the superconducting phase 0(r) passes through 
TT twice, separated by a time interval of T5. Therefore, 
the zero energy eigenvalues at (/) = tt split to finite ener- 
gies SX = ±EMe-^^^'/2. When El > 2Ej/{37r^) and 
hence ^5 < tt, the imaginary time interval of T5 can be 
readily evaluated from 



V^j(l - cos (/)) + El^^ - A^^El ' 

(33) 

For El/Ej <C 1, we can ignore the contributions from 
El/Ej from the integrand. Thus, this integral can be 
approximated by 




2^2EcEjn - I 

J 6 



I , , , ,1 = — In tan ^ 
|sin(^/2)| 4 



(34) 



with 4>b « 2t:\J2El/ Ej. By droping the constant terms, 
we have 

1 



HEj/El). 



(35) 



In Fig. [To) we show the numerically evaluated T5 as a 
function of El/Ej with Ec = Ej = 1. The approxi- 
mated T5 in Eq. ( [35| ) is in good agreement with numerical 
results. 

From Fig. [sjc), we observe that most eigenvalues re- 
main unchanged in the presence of an instanton despite 
the appearance of zero modes. As the zero energy modes 
split to 



SX = ±:^(^^/^^)£M/(2y2£ZJB7)^ 



(36) 



the tunneling rate is changed by the ratio of the deter- 
minant of the fermionic kernel in the presence and in the 
absence of the bounce. This ratio is dominated by 



F (X 



\sx\ 



y^det[Lffi] Em/2 \E, 



Er 



Em/{2V2EcEj) 



(37) 
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This result is in perfect agreement of the suppression 



of relaxation rate e [given in Eq. (30)] due to the 



presence of fermions. 

Let us now discuss the relevant energy scales and 
the experimental feasibility of such a system. First, 
for the bouncing event to cross phase tt, it requires 
El < 2Ej/(37r2) - OMlbEj. We also need El > 
£^M/(47r^) to make phase slips of 27r energetically possi- 
ble. Therefore, we require a system satisfying the condi- 
tion > > Em/ {4^71^)- Finally, we need Em ^ Ec 
to make the dependence on El observable as it requires 
that the exponent in Eq. (37) is of order unity. 



We shall seek an experimental construction with a 
large Ej/Ec ratio such that the energy scale hierarchy 
can be realized. In general, a Josephson junction with 
Ej ^ Ec can be made out of a Nb/AlOx/Nb junc- 
tion. A typical critical current density of such a junc- 
tion with insulating layer thickness 1 ~ 10 nm is in the 
range of jc = 10 ~ 1000 A/cm^, see Ref. 40, For a 
junction of area 10~^ cm^ with critical current density 
jc = 20 A/cm^, we can estimate the Josephson energy 
by Ej = ^o^/(27r) ?^ 5 K. With the thickness of the 
insulator at 5 nm, the expected capacitance of such a 
junction is about 18 fF and leads to a charging energy at 
Ec ^ 200 mK. For a semiconductor wire in contact with 
Niobium, Em can be of the order of 0.1 — 1 K as the super- 
conducting critical temperature Tc ~ 9.2 K for Niobium. 
Here, we will assume that Em ^ 0.5 K, which gives the 
exponent in Eq. ^ as E m/ (2^j2E^Wj) - 0.18. Fi- 
nally, we need a relatively large inductance L > 12 nH to 
satisfy El < 0.067bEj. Such values of inductance can be 
achieved with a larger ring or with a more complicated 
desie; In the followine;, we will show that the same 
physics can be accessed in a much simpler setup without 
inductance at all. 



B. Current biased geometry 

The second geometry we consider is that of a Josephson 
junction on a topological superconducting wire, and we 
pass a supercurrent Is through the wire. The effective 
action is then 



Jo 



Sefi = I dr 
Em 



1 1 



2 8Ec 



(t)2 + ^j(1-cos0(t)) (38) 



±^cos(0(t)/2). 



h . 

Ye' 



(39) 



which has a titled doubly periodic washboard potential. 

In the case where there is no supercurrent applied, 
/g = 0, the system relaxes to a stationary state where 
the superconducting phase difference is pinned to a mul- 
tiple of 47r. Successively, driving the system with an ex- 
ternal current of size Ig tilts the potential. The system 
is trapped in a metastable state having the possibility to 
tunnel through the potential barrier out of the local min- 
imum. Employing the same analysis as in the previous 



subsection, one can show that the presence of the exter- 
nal current plays a similar role as the inductance term 
in the ring in particular it makes the potential minima 
separated by 27r tilted, thus giving the system an incen- 
tive to tunnel and thus lower its energy. In particular, 
the effect of the bias supercurrent on the relaxation rate 
r is given by (37) with 



El ^ 



4e7r 47r^ 



(40) 



Differently from the previous setup, the phase of the 
current-bias wire after tunneling enters a so-called run- 
ning state which means that the wir e tur ns resistive, es- 
sentially switching to a normal stat d^^^^M After turning 
off the bias current, the superconducting phase retraps 
in one of the minima due to dissipation given by a small 
shunt resistor. 

The experimental determination of the relaxation rate 
r thus goes along the following line. First, the current- 
bias is turned off and the wire is prepared in its ground 
state. Then, the bias is turned on to a value Is on a 
timescale Ton <C F"^. The time difference between the 
event of turning on the current bias and the switching 
of the voltage to a finite value is a direct measure of the 
inverse phase-slip rate F"^. Repeating the experiments 
for different values of Is the predicted power law (37) 



can be tested and thus the suppression of the quantum 
phase-slip rate due to the zero-mode when lowering Is in 
the regime Air'^Ej Is^o ^ Em can be confirmed. 

Thus far, we have assumed that the initial state (be- 
fore we turn on the bias supercurrent) corresponds to 
the Josephson junction localized in the deeper well of 
the doubly periodic potential. Alternatively, we could 
prepare the Josephson junction so that it is localized in a 
random well (e.g. by driving it). With this type of initial 
condition, there will be two relaxation rates, correspond- 
ing to the two types of wells in the doubly periodic po- 
tential. Thus, the experimentally observed distribution 
of waiting times should be bimodal. 



VI. CONCLUDING REMARKS 

We investigate phase slips in topological superconduct- 
ing wires. Unlike in conventional superconducting wires, 
phase slips in topological superconducting wires occur 
in multiples of 47r as opposed to multiples of 27r. Our 
original motivation for looking into this problem was to 
understand the effects of phase-slips in topologically pro- 
tected qubits made up of conventional and topological 
superconducting wires. As phase-slips are non-local per- 
turbations, they can cause decoherence of a topologically 
protected qubit. 

The fact that phase-slips in topological wires occur in 
multiples of An is well known. Indeed, by integrating out 
the fermions, one finds that the effective action for the 
phase is 47r periodic. We show an alternative explanation 
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of this fact by a beautiful analogy to spontaneous sym- 
metry breaking of the theta vacuum in quantum chro- 
modynamics. For the case of QCD, t'Hooft found that 
in the background of the instanton of the gauge field, 
there is a zero mode in the fermionic determinant P^^^ 
This zero mode results in the vanishing of the transition 
rate between configurations of the vacuum with different 
winding numbers. Similarly, we find that in the back- 
ground of a 27r phase slip, the fermion determinant con- 
tains a "hidden" zero mode, that results in the vanishing 
(suppression) of the rate of 27r phase slips. 

Returning to the question of decoherence of qubit de- 
vices, we show that phase-slips in the topological wires do 
not cause decoherence, as they occur in multiples of 47r. 
However, the qubits are susceptible to decoherence from 
27r phase slips in the conventional superconducting wire 
segments. Phase-slips near the junction of conventional 
and topological superconducting segments can also occur 
in multiples of 27r, and likewise result in decoherence. 

To investigate the 47r periodic nature of topological 
superconductors we propose two types of experiments. 
First, we suggest looking for changes of magnetic flux 
through a ring made up of a topological superconduct- 
ing wire by multiples of two flux quanta as opposed to 
one flux quanta for conventional superconducting rings. 
Second, for a current biased wire, we suggest looking for 
voltage spikes corresponding to 47r phase slips as opposed 
to 27r phase slips. 
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Appendix A: Weak link model 

In this appendix we derive the effective fermionic ac- 
tion for the weak link geometry. We start from the Kitaev 
model ([T]), and extend it by letting the magnitude of both 
ti^i-^i and A^^^+i to vary from site to site 

N-l 

H = -Y1 (^i,i+i4+iCi + h.c}j 

N-l 

J2 (Am+i4+ic1 + h.c.) . (Al) 



A=l A=s e'* A=e'* 



i=l 



Next, we set A^^^+i = ^^^^+16*'^^'^+^ on all links. To model 
the weak link geometry, we need to (1) make a weak link 



original t y — " — V 

fermions 



Majorana 



fermions 



link 

fermions 



^1 ^2 ^j+l 

A A AAA A A 

000 0'"'0 0000 o»— o 000 



Li Ri L, 






(^^^^^^ 0000 (^^^ (^^^^^^ 0000 (^^^^^ 
fl fj fj+1 fN-1 



FIG. 11. Schematic of the Kitaev model extended to con- 
tain a weak link. The top row, labeled "original fermions" 
shows the original model and indicates the values of the pair- 
ing function on the links. The weak link, going from j to 
j + 1, is highlighted in blue. The middle row, labeled "Majo- 
rana fermions" , indicates the transformation to the Majorana 
fermion basis. The bottom row, labeled "link fermions" shows 
the final transformation to the link fermion basis. The link 
fermion associated with the weak link is highlighted with pur- 
ple shading. 



and (2) set the phase of the pairing field to be different 
to the left and to the right of the weak link. Therefore, 
we (1) set ti^i-^i = 1 on all links except the weak link, on 
which we set tjj-^i = s with s < 1; (2) we set the phase 
= for z < j and (/>i,i+i = (/> for z > j, see top row 



of Fie;. 11 



Having defined the model, we can obtain the eigenspec- 
trum following the schematic steps illustrated in Fig. [TT 
First, we rewrite the Hamiltonian in slightly more con- 
venient form 



N-l 



+ 1^ 



. (A2) 



Now we introduce "right" {R = R^) and "left" {L = L^) 
Majorana operators at site i as 



equivalent to 



4 = ^ [Ri^iH 



g+j</)i,,+i/2^t _|_ g-i0i,j-|.i/2g 



(A3) 



(A4) 



Ri = 

which allow us to write 
H = ^ iRiLi-^i 

+ isRj [sin(^/2)i?_^-+i + cos(0/2)L_^-+i] 



Note at this point that the two Major anas at the ends of 
the wire, Li and i?Ar, do not appear in the Hamiltonian - 



13 



they do not couple to anything, and therefore constitute 
the two Majorana zero modes at the two ends of the 
wire. We can recombine the Majorana operators to form 
fermion operators on the hnks 



for which 



=-[l + ii?,L,+i] 



where we used the fact that i?| = Lf 
terms of the Hnk fermions we get 



{ci,4} 



(A5) 



(A6) 



1. In 



sW) (2/j/^. - l) 
-z.sin(0/2)(/,l+/,) (4i+/,+i) . 



5 COS 



The two Majoranas at the ends of the wire, Li and 
i^AT, may be combined into a single complex fermion 
/o = ^ [Rn — ^^i], which does not appear in the Hamil- 
tonian. After transforming to the link operators, we find 
a Hamiltonian that is almost diagonal. The exception 
being the terms involving fj and /j+i operators in the 
vicinity of the weak link. Explicitly, the non-diagonal 
part of the Hamiltonian is 



scos(0/2) issm{(j)/2)/2 is sin(0/2)/2 \ 

-5cos(0/2) issm{(j)/2)/2 i5sin((/)/2)/2 | 

-issm{(j)/2)/2 -is^m((j)/2)/2 1 

-issin(^/2)/2 -issin(^/2)/2 -1 / 



/ 

V /j+i J 



(A7) 



The non-diagonal part of the Hamiltonian may be readily 
diagonalized via a Bogoliubov transformation with eigen- 
values 



252cos(0)). (A8) 



Here, we denote two positive eigenvalues (at ^ = 0) to 
correspond to the annihilation operators and and 
the two negative eigenvalues to the creation operators cj^ 
and djy. Away from (/) = 0, we label the operators such 
that their eigenvalues evolve smoothly. Putting all these 
considerations together, the Hamiltonian becomes 



H 



E 

=l..i-lj+2..iV-l 



(2///, - 1) 



(A9) 



Ox 



(2/0V0 - 1) 




2s^ cos(^) 



(24 



c - 1 



2s^ cos((/)) 



{2Syj 



1 



where the first line corresponds to all the link fermions 
except fj and /j+i, the second line corresponds to the 
zero mode on the ends of the wire, and the final two lines 
represent the link fermions, fj and /j+i, around the weak 
link in their diagonalized basis. 



To better understand the fermions on the links j, j + 1 
and j + l,j + 2, we plot the eigenvalues as a function 
of the phase difference across the weak link in Fig. [T2j 
The figure shows two distinct types of eigenvalues: (1) 
The eigenvalues correspond to are near zero energy 
with eigenfunction localized mostly on the link j, j + 1. 
(2) The eigenvalues correspond to duj are near ±1 with 
eigenfunction localized mostly on the link j + 1, j + 2. 
Case (1) corresponds to the eigenvalues that cross zero 
at TT and Stt. The crossing indicates that the creation and 
annihilation operators and cj^ switch after the phase 
rotates by 27r. As we make the weak link stronger (by 
increasing s), the two types of eigenvalues approach each 
other at ^ = 0,27r,47r. In the limit 5 = 1, we return to 
the case of a continuous wire, and the eigenvalues touch 
at multiples of 27r. 



Using our understanding of the adiabatic evolution of 
the eigenvalues, we can choose the correct branch of the 
square root in writing the penultimate term of the Hamil- 
tonian (|A9[) (i.e. the solid line in Fig. 12). Thus, expand- 



ing to second order in 5, we can write the Hamiltonian 
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(a) s = 0.25 



(b) s = 0.75 




FIG. 12. Spectrum of the four eigenvalues of (A7) as a function of the phase difference (f) across the weak link for three values 



of s parameter as indicated. The traces correspond to adiabatic evolution of the eigenvalues, two eigenvalues indeed cross at 
TV and Stt. Thus we see that the spectrum has both a 27v periodic component (the two eigenvalues at e ~ ±1) as well as a Att 
periodic component (the two eigenvalues at e ^ 0). 



as 



i=l..j-l,j+2..N-l 

+ 0x (2/0V0-1) 

+ scos(^/2) (24c^ - 1) 



(AlO) 



1 + — [1 - cost 



1 



What does the ground state of this Hamiltonian look 
like? Starting from the state |0) with no q fermions we 
obtain two degenerate ground states 

\GSi) = foc^d^l[m (All) 

i 

\GS2) = c^d^Ylm (A12) 

i 

corresponding to the zero mode at the ends of the wire 
filled and empty. 

Thus far, we have implicitly assumed that our opera- 
tors are phase dependent, i.e. = c^(^). To under- 
stand the evolution of the ground state, we can relate 
the phase evolved operators to the operators at = 0. 
Under the rotation of the phase (j): ^ 27r the operators 
(as defined at ^ = 0) evolve as: 

/.(«)- {4!?i.;5i, (A13) 

/o(0)^-/ot(0) = /o(27r) (A14) 
c^(0) ^4(0) =c^(27r) (A15) 
d»(0) ^ -d^iO) (A16) 

Using these rules we can easily work out the evolution of 
the ground state: 

|G5i) ^ (-l)^-Vo^(0)c|„(0)d»(0)n/i(0)|0) , (A17) 

i 

\GS2) ^ i-ir-^-'cliO)dUO)l[fm\0) ■ (A18) 
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FIG. 13. Generic fermion spectrum as a function of the phase 
slip angle in a topological wire showing three types of modes: 

(1) zero mode associated with the Majorana fermions at the 
ends of the wire (red dotted line along the horizontal axis), 

(2) Atv periodic mode associated with the fermion localized on 
the phase slip, and (3) a set of 27r periodic modes associated 
with the un-localized fermions. 



Applying the operator transformation rules twice, we find 
that a 47r phase slip returns the system to the original 
configuration. 

Thus far, we have obtained the effective Hamiltonian 
for the very special case of the Kitaev model at |A| = t. 
However, we point out that the form of the Hamiltonian, 
and of its spectrum, is generic: containing three ingredi- 
ents (1) a zero mode associated with Majoranas at the 
ends of the wire (2) a single mode that corresponds to a 
Fermion localized on the weak link that is 47r periodic, 
and (3) a set of 27r periodic modes. In particular, this 
form of the spectrum is not sensitive to adding small dis- 
order, the relaxation of the condition |A| = t, nor the 
spreading out of the "weak link" over several links of the 
model. We plot the spectrum for several generic cases in 



Fig. 13 



To understand the evolution, we need to obtain the 
effective action. The imaginary time (Euclidean) action 
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for the Kitaev model is 



S 



f 

Jo 



dr 



N 



Y.c\{T)drC,{T)-H 



(A19) 



We now focus on the first term, and perform on it the 
same transformation that we used to diagonalize the 
Hamiltonian. Transforming first to the Majorana fields, 
we have 



[Ri + iLi] dr- 



[Ri - iLi] 



Ridq-Ri -\- Lidq-Li -\- iLidq-Ri — iRidq-Li 



- idr(t){l + iLiRi) 



(A20) 



where we have allowed for the possibility that the su- 
perconducting phase (j) changes with (imaginary) time. 
Using the anticommutation of Majorana fields, and inte- 
gration by parts we find 

LdrR - RdrL = LdrR + {drR)L = {L, drR} = , 

(A21) 

from which we get to 



cKr)drC^{r) 



RidrRi + LidrLi + i9r^(l + iLiRi) . 

(A22) 



Next we want to change fields to the link fermions. Using 
(A21) again, we can show that 



//(TR/,(r) 
which yields 



1 



[RidrRi + Li+idrLi+i] , (A23) 



cl{T)drC,{T) = f}{r)drfi {r)+\idr4>{i + iLiRi) . (A24) 



Now we also need 



iLiRi 



/i+l fi + // 



where we have used 



(A25) 



(A26) 



We see that allowing the phase to change 
dynamically introduces an off-diagonal term 



/i+l /i+l fi + fi 



in the link fermion 



action that corresponds to exciting pairs of link fermions 
on the neighboring sites. If is much smaller than 
the system gap, this process will be forbidden, and this 
term can be safely ignored. We shall assume that we are 
indeed working in this limit (e.g. due to extrinsic factors 
like the capacitance that dictates how fast the phase can 



evolve). Although the off-diagonal term has no effect on 
link fermions away from the weak link, what about the 
fermions and d^? Applying the transformation used 
to diagonalize Eq. (A7), we find that we must again 
excite pairs of quasi-particles (i.e. terms of the form 
cjy(ijy), which is again forbidden by energy conservation. 

Combining these considerations, and integrating over 
all link fermions except the fermion, we arrive at the 
effective action for the weak link 



^eff 



icos(^/2)(2cT, 

o2 



1) 



(A27) 



1 



[1 



Appendix B: Discretization of fermionic kernel 

In this Appendix, we establish an alternative dis- 
cretization scheme that avoids the second fermion dou- 
bling for the Lagrangian density 



dr^ A(r) 
~ ^ dr- A(r) 



(Bl) 



in Eq. ([17|), where A(r) 
ferential equation. 



^MCOs^(r). Hence, the dif- 



(lLS, + a,A(r)) 



u{r) \ _ 
v{r) 



A 



u{r) 
v{r) 



(B2) 



can be solved numerically. Then, the determinant of the 
kernel is simply given by the product of all eigenvalues, 
A, in the range of r = (0, /3) with /3 ^ oo. 

The fermionic kernel appearing in the eigenvalue equa- 
tion, Eq. (B2), can be brought into a hermitian form 



by multiplying from the left with —icjy^ namely Hf = 
—icjyCf. This results in the eigenvalue problem 



{-iCFydr + (7x/^{t)) 



v'{t) 



w / u'{t) 



(B3) 



Such differential equation, in turn, describes the contin- 
uous theory of the polyacetylene with identification of r 
as the one-dimensional coordinate and A(r) as the mass 
field. As det[icr2y] = 1, we have det[>C/] = det[i^/] = 
HA'. 

It is well known that a zero energy mode in Hf will ap- 
pear and localize at a soliton where the mass field changes 
sign. In our system, the mass field changes its sign when 
the superconductor phase alters from (/) > tt to ^ < tt 
or vice versa. Therefore, we expect a vanishing determi- 
nant, as one of A' is zero, when such phase change occurs 
in the imaginary time range r = (0,/3). The appearance 
of such zero energy modes, however, becomes illusive in 
Eq. (B2). Hence, a proper discretization scheme to make 
it a difference equation becomes handy for understanding 
the effect of those zero modes. 

Let us first review how to obtain the continuum limit 



of Polyacetylene, Eq. (B3), from the lattice model. Then, 
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the discretization scheme for Eq. (B2) can be inferred by 
reverse engineering the corresponding parameters. Here, 
we win fohow closely the derivation in Ref.[43|and assume 
open boundary conditions. 

We consider the following one-dimensional lattice 
Hamiltonian on 27V sites 



j=2n-l 



j=2n 



where aj{aj^) and bj{bj^) are annihilation(creation) op- 
erators at the odd and even lattice sites, respectively, 
and n > 1 are integers. Here, tj-^ij are real hopping 
amplitudes between site j and j + 1 and take the form 



t 



^0 - 7(^. 



(B5) 



where yj is the displacement of the atoms at site j and 
7 is a constant governing the variation of the hopping 
strength due to the displacement. The lattice constant 
is ^/2 and hence the distance between two adjacent a (or 
b) type atoms is S. As the lattice displacements reflects 
the induced Peierls instability, they take the form yj = 
{ — lyrjj. Then, the hopping amplitudes become 



t 



(B6) 



Here, r]j can be a function of space, but if rjj are uniform 
the hopping amplitudes are simply alternating between 
even and odd sites. 



To match the continuum theory with Eq. (B3), it is 



convenient to introduce the gauge transformation 

a,^(-l)^(i)e-*-/%, b,^{-iy/'e^-/%. (B7) 
and then the lattice Hamiltonian becomes 



j=2n-l 



j=2n 



We note that Eq. (|B4|) and Eq. (|B8|) differ only by the 
gauge transformation and have exactly the same spec- 
trum. 

By arranging operators with 2N lattice sites in a vector 
form 

^ = (ai, as, . . . , a2N-iMi ^4, • • • , ^2iv)^ , (B9) 
the lattice Hamiltonian becomes 



Q 

3^ 



(BIO) 



(a)i 



10"^ 
10-^ 
10-^ 
10-9 



cosh[/9/2]-2 

|L/,27r|/|^/,o| 

— zero modes 
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FIG. 14. Spectrum using the modified discretization scheme 
Eq. (B19), that avoids the additional fermion doubling, (a) 
Ratio of determinant, same as Fig. |6] using modified dis- 
cretization scheme, (b) and (c) same as (b) and (c) of Fig. [s] 
using modified discretization scheme. 



where Q is an A/" x A/" matrix and has the form 

\ 

io -to U U ■•• 

V ; 



-^0 











to 


-to 











to 


-to 











to 


-to 



I 

/r]2^Vi 

^5 + ^4 ^6 + ^5 

V ! ! : ••. 



(Bll) 



In what follows, we will show that the first sector of Q 
corresponds to kinetic energy while the second sector of 
Q corresponds to the mass field. 



To derive the continuum theory of Hamiltonian (B8), 
we first identify the operators with continuums fields 



(B12) 



a, = {sY/'u(^{j + i)iyb, = {sY/'v^j 

Inserting this into Eq. ( |B8| ), we have 

Hp = + to5j2vHri5]{U[{n+l)S]-U[n5]) 

n 

- to5j2uHn5]{V[n5] - V[{n - 1)5]) 

n 

+ iS Y,(V2n + mn-i)V''[nd]U[nd] + h.c. 

n 

+ 15 Y,(V2n+i + V2n)V^ [n5]U[{n + 1)5] + h.c. 
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If we assume that all [/, V, and rj are smooth functions 
in x-direction and expand the Hamiltonian to the linear 
order in J, the Hamiltonian becomes 



Hp=^ toS^ V^{nS)U\nS) - U\nS)V\nS) 

n 

+ 475 ^ 7?(n<5)0(n<5)/7(n<5) + /i.e., 



(B13) 



We notice that the discretization of the mass field is 
spread out over several matrix elements. As Eq. (B2) 
and Eq. (B3) differ by a multiplication of —icy^ the dif- 



ference operator of kernel C is given by 



where 



U\nS) ={U[{n + 1)S] - U[nS])/S 
V\nS) ={V[nS] - V[{n - 1)S])/S. 



(B14) 



By taking the continuum limit, we replace 
fdx, U{n5) U{x), V{n5) V{x), V\nS) d^U{x), 
V'{n5) dxV{x) and r]{n5) Tjix)^ and obtain 



+ 47 / dxri{x){V\x)U{x) + U\x)V{x)). 



(B15) 



The continuum Hamiltonian can be rewritten with the 
vector operators x(^) = {U{x)^V{x))^ as 



Lf {icFy (g) tN)l-LL = 





-C 



(B19) 



dxx\x)[{toS){-id)ay + 4-fr]{x)a:,]x{x). (B16) 



(B17) 



This matches the form of Eq. (|B3|) with 

toS = l 
4777(0;) = A{x) 



to = 1 A 

777(0;) = A(o:)/4. 



By defining 777^- = A(j|)/4 = Aj/4, we have the 
lattice identification of the continuum operator H = 
—idxCTy + A{x)(Jx as in Eq. (BIO) and 



/-I •••\ 
1 -1 ••• 
1 -1 ••• 

/ A2 + Ai 
A3 + A2 A4 + A3 
A5 + A4 A6 + A5 



J 



(B18) 



where In is a N x N identity matrix. 



We compare the results for fermions in the background 
of a phase-slip followed by an anti-phase slip using the 
original discretization scheme of Sec. IV and the alter- 
native discretization scheme discussed in this appendix. 
Using the alternative discretization scheme, we find that 
in the continuum limit the value of the fermion determi- 



nant after taking a square root (Fig. 14 'a)) is identical 
to the one obtained in the original scheme after taking 
the fourth root (Fig. |6|. However, the spectrum of the 
untransformed Lagrangian (Fig.[l4jb)) does not have the 
nice properties found in the original scheme (Fig. [5|b)) 
due to the explicit mixing of the left and right moving 
components. The number of hidden zero modes indeed 
decreases from four (Fig. [sjc)) in the original scheme to 
two (Fig. 14 'c)) in the alternative discretization scheme. 
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